Distribution of Brownian particles in turbulence 
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We consider Brownian particles immersed in the fluid which flow is turbulent. We study the 
limit where the particles' inertia is weak and their velocity relaxes fast to the velocity of the flow. 
The trajectories of the particles in this case have a strange attractor in the physical space, if the 
particles' diffusion is neglected. Under the latter condition the singular density of the particles was 
recently described completely. The analysis was done for real turbulence and did not involve the 
flow modeling. Here we take the diffusion into account showing how it modifles the statistics. The 
analysis is performed also for real turbulence. Experimentally testable predictions are made. 
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Lately the problem of spatial distribution of inertial 
particles in turbulent flows enjoyed much attention of 
the researchers pl-fTsj. The subject became timely after 
the breakthrough in the theoretical understanding of the 
Lagrangian motion of particles in the flow occurred [l^ . 
While the more traditional subject of tracers that follow 
the flow is quite complete by now [l^ . the understand- 
ing of the behavior of inertial particles is still insufficient. 
Since all real particles have inertia, and the effect of iner- 
tia was shown to be a non-perturbative phenomenon that 
often cannot simply be neglected, then the study of the 
subject is a very important practical question. Theoret- 
ical advancement was made mainly for the case of small 
Stokes number, where the inertia is weak and the parti- 
cles "almost" follow the flow. Even in this limit of small 
Stokes numbers, the particles' distribution is highly non- 
trivial. Particles' deviations from the surrounding flow 
accumulate with time, bringing particles to a strange at- 
tractor in space. This attractor is multi-fractal and its 
statistics was recently described completely in [l^, see 
also jl6i] . The description was made for the real turbu- 
lent flow and it did not involve any modeling. Within this 
model-independent approach, the statistics of turbulence 
is considered as given but not known. The final answers 
are expressed with the help of statistical characteristics of 
turbulence. To apply such kind of results one should first 
estimate or measure the corresponding statistical prop- 
erties of turbulence, and then substitute them in the an- 
swer for the particles. This approach was first used in 
[3| where the correlation codimension of the particles' 
attractor was expressed in terms of different time corre- 
lation function of the Laplacian of turbulent pressure. It 
was continued in [l^ where all the rest of quantities, be- 
sides the correlation codimension were obtained. Since 
the statistics of turbulence is largely unknown 17[ , then 



obtaining such results demands understanding of the uni- 
versal features of the particles' behavior in the flow, that 
are independent of the details of the flow statistics. 

In this Letter we continue the model-independent ap- 
proach outlined above to include the Brownian motion of 
the particles. We again consider the case of small Stokes 
numbers, adding the Brownian motion to the underlying 
equations of the particles' motion. When the Brownian 



motion, or diffusion, is taken into account, the density 
of the particles in the steady state becomes an ordinary 
smooth, random field. We describe the statistics of that 
field and make a number of predictions that are testable 
experimentally. 

The idea that particles' inertia leads to inhomogeneous 
spatial distribution dates back to the seminal paper by 
Maxey 3- It was observed that, due to the inertia, heavy 
particles are pushed out of the vortices and hence they 
will not distribute uniformly in the flow, like the inertia- 
less particles. 

Our analysis is based on the recent finding of a uni- 
versal description for the singular steady state density 
of the weakly compressible dynamical systems [l^ . The 
particles' motion, though governed by Newton's law, ad- 
mits an effective description in terms of a velocity field 
in space, where inertia is described by a small compress- 
ible correction to the incompressible velocity of the back- 
ground turbulent flow. This correction leads to a small 
disbalance of trajectories going in and out of space re- 
gions, which accumulates over a long time to a big effect. 
Thus, compressibility is a singular perturbation. The 
adequate treatment was performed in [l6| . It was shown 
there that at sufficiently small scales the inertia is always 
important, however small it is. When diffusion is taken 
into account, the steady state density is a smooth func- 
tion. This function is significantly different from a con- 
stant only if the diffusion coefficient of the particles is not 
too large. Hence we will assume that the Prandtl num- 
ber is large, that is the fluid viscosity is much larger than 
the diffusivity of the particles. It will be clear that large 
Prandtl numbers are necessary if any significant fiuctu- 
ations of the density are to be present in the considered 
case of small Stokes numbers. 

We pass to the analysis. We assume that the density 
of the particles is not too large, so the hydrodynamic in- 
teraction between the particles and their back reaction 
on the flow is negligible. Then the analysis is effectively 
the one of the single-particle motion. We consider a small 
spherical particle with the radius a and the material den- 
sity pp suspended in the fluid with the density p and the 
kinematic viscosity v. The fluid flow u(t, x) is assumed 
to be incompressible. We first consider the equations 
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that neglect the Brownian motion of the particles. The 
Newton law governing the evolution of the particle's po- 
sition x{t) and the particle's velocity v{t) is assumed to 
have the form 



dv 



V — u[t, x{t)] 



(1) 



where r = (2/9)(pj,/p)(a^/j/) is the Stokes time. Thus 
we assume that all the forces besides the drag can be 
neglected 0, • Above the Kolmogorov scale 77 is 
assumed to obey 77 a. With no loss we set the total 
volume and the mass equal to one, so the spatial average 
of the particles' density n obeys (n) = 1. After transients 



v{t) 





■ t-f 


/ exp 




/ —00 


T 



dt' 

u[t\x{t')] — . (2) 
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We assume r is much smaller than the smallest time-scale 
of turbulence, which is the viscous time-scale t,,, so the 
Stokes number St = r/i^ <C 1. Then we can substitute 
the derivative in the integrand by its value aX t' = t 
so v(t) u - T [dtu + {u ■ V)u]. Thus at S'i < 1 the 
particle's velocity is determined uniquely by its position 
x{t) in space and one can introduce the particle's velocity 
field v{t,x), see 

x{t) = V [t, x{t)] , V = u — T [dfU + {u ■ V)m] . (3) 

In the zero inertia limit St — > the particles follow the 
incompressible mixing flow of turbulence x = u [t,x(t)] 
and in the steady state they are uniformly distributed in 
space, so their steady state density equals one. This 
behavior is characteristic of small dye particles. However, 
at a small but finite St, the small correction v ~ u gives 
the particles' velocity field a finite compressibility [2] 



-rV • (m • V)m 7^ 0. 



(4) 



The physical significance of the non-vanishing diver- 
gence can be understood by considering the inertial 
outward motion of heavy particles in a centrifuge, see 
[11) S 0, [ill ) ^■nd also [l^ for the relation quantifying the 
anti-correlation between the vorticity and the particles' 
density. The steady state density Ug of the particles mov- 
ing according to Eq. ^ is singular. It is supported on the 
strange attractor - a multi-fractal set in space which is 
asymptoted by the particles' trajectories at large times, 
see [i^l for the complete description of that set. The par- 
ticles' density has log-normal statistics and it is described 
completely by the pair-correlation function 
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(5) 



where angular brackets stand for spatial averaging, and 
< Cky <C 1 is the Kaplan- Yorke dimension of the 
attractor [l^. The latter is estimated as St'^ times a 
positive power of the Reynolds number [3, IB [Hi • The 
divergence of (n^) corresponds to the singularity of the 



density described above. Note that by Cky ^ 1 there 
are no significant fiuctuations of the density already 
at scales x <^ i], which is inferred by the saturation 
{ns{0)ns{x)) « 1 = (us)-^ that holds at a: ^ 77. The 
presence of the separation scale / <^ rj where there are no 
significant fiuctuations of the density is the most funda- 
mental property of the weakly compressible regime, cf. 

The above results however cannot be applied directly 
to the real particles for which the Brownian motion is 
present always. Here we include the Brownian motion 
into account, studying how it changes the steady state 
density. Our starting point is the equation of motion 



dx 
'dt 
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where k is the diffusion coefficient. Below we assume that 
the Prandtl number Pr = v/k, is much larger than one, 
Pr » 1. 

It is clear that the diffusion will smoothen the sin- 
gular steady state density described above by making 
the particles diffuse from the spatial regions, where they 
would be localized otherwise by the combined action of 
the fiow and the inertia. The particles' density n obeys 
the Fokker-Planck equation 



le pa 

m 



dtn + W ■ [vn] = kV' 



n, 



(7) 



where the diffusion coefficient stands near the higher- 
derivative term, suppressing non-smooth realizations of 
n. Further, we note that because diffusion only smears 
the inhomogeneities making the density smoother, then 
it follows that there are no significant fluctuations of the 
density a,t I rj also with diffusion. This property is at 
the center of the consideration below. 

The key tool of our analysis is the Green's function 
G{x\t,x') that expresses the density n{0,x) in terms of 
the density n{t,x) at t < according to 



n(0, x) 



J G{x\t,x')n{t,x')dx', G{x\0,x') 
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We analyze G{x\t,x') as a function of its initial argu- 
ments t and x' . As a function of t and x' it satisfles the 
equation 



dtG + v{t, x') ■ V^,G = ~K,Vl,G, 



(8) 



which is the Hermitian conjugate to Eq. ([7]), see [l^. We 
follow the analysis of 3]. At k = the Green's function 
is supported on the particle's trajectory q{t, x) obeying 



dt 



^ v[t,q{t,x)], q{0,x)^x. (9) 



One has 

G{x\t, x') = cxp w I^') <?(*') x)]J ^ - 9(^) ^)1 • 



3 



Due to the diffusion the particles that are not located at 
q{t,x) at time t can also arrive at x at time 0. This is 
summarized by the diffusive smearing of the support of 
G{x\t,x') near x' = q{t,x). If we introduce G{x\t,x') = 

exp ^- /° w [t', q{t', x)]j ^'[i, x' - q(t, x)] then obeys 

the equation that follows from Eq. and reads 

dt^^ + iv [t,x' + q{t,x)] - v[t,q{t,x)]) • V^-* 
^^KVl,'^-w[t,q{t,x)]^ (10) 

We note that at not too large t, that will be deter- 
mined below, the dimensions of the support of ^[i, a;] 
are much smaller than rj, so that one can approximate 
v[t,x' + q{t,x)] — v[t,q{t,x)] sa a{t)x' where aij{t) = 
WjVi[t,q{t,x)] and tra — w[t,q{t,x)]. Making the 
Fourier transform in x' we obtain 

dt'i' - cTtjh-Vj^ = nk^^, *(t = 0) = l. (11) 

The solution to the above equation can be written as 



*(t, fc) = exp 



-kilij{t)kj 



where the symmetric matrix I{t) obeys 
dl 



dt 



al + la' - 2k, I(t = 0) = 0. 



(12) 



(13) 



The matrix I{t) admits an explicit expression in terms of 
the Jacobi matrix Wij{t, x) = Vjqi{t, x) of the backward 
in time flow q{t, x). This matrix obeys W = aW, so that 



I{t)=2KW{t) j W-'^{t',x)W-'^'\t\x)dt'W\t). 
It follows that 

G{x\t,x') = exp(^~ w[t',q{t',x) 



dk 



exp 



ik ■ {x' - q[t,x]) - -kl{t)k 



,(14) 



where d is the dimension of space (in the physically rel- 
evant situation d = 2 or ci = 3). The above answer 
coincides with the one of Q with a slightly different no- 
tation. The equation makes it clear that I(t) is the iner- 
tia tensor of the cloud of particles localized at t = at 
a; = 0, cf. [20]. Due to the diffusion the cloud's particles 
spread around q{t, x) over an ellipsoid which axes are de- 
termined by the eigenvalues of y^I(t). The largest axis L 
of this ellipsoid grows as L{t) = Idif exp[\Xdt\] where we 
defined the diffusive scale Idif as Idif = y/Kt^. This scale 
arises as follows. The cloud of particles "initially" local- 
ized at X starts spreading backward in time. The spread 
is initially completely diffusive as the particles moving 
according to Eq. ^ have initially very close velocities. 
In particular, in Eq. (jl3p at small t, the diffusion term 
2k is the largest and Iij{t) « 2K\t\Sij. However at times 



t of order the diffusive term becomes of the order of 
a I terms. The subsequent growth was described in j20| . 
The application to the weakly compressible case consid- 
ered here gives that at < ^ the scale L{t) grows due 
to the chaotic exponential separation of the trajectories 
backward in time. This separation occurs at the princi- 
pal exponent of the flow backward in time that in our 
case coincides with the absolute value of the d— th Lya- 
punov exponent [2l[ of turbulence Xd ■ This is because the 
flow is weakly compressible and to the leading order the 
Lyapunov exponents of v coincide with those of u. The 
latter flow is incompressible, so the principal exponent of 
the flow backward in time is |At;|. 

We now can determine the domain of applicability 
of Eq. p4p . This equation was derived assuming the 
linear size of the support of ^ is much smaller than 
T]. This means L{t) ^ 77, i. e. Eq. (|14l) applies as 
\t\ « \Xd\~'\n{v/ld^f), cf. [3j] 

We are now ready to find all the single point moments 
of the density. We use the presence of the separation scale 
I ^ rj at which there are no appreciable fluctuations of 
the density. Consider 

n(0,0) = exp w[t',q{t',x)] 

X J 'i/[U,x' - q{U,x)])n{U,x')dx' (15) 

where = ~\Xd\~^ ^nil/ldif)- The function in the 
integrand (which has a unit space integral) performs an 
effective averaging of n(i*,a;) over a region with linear 
size I. Since there are no fluctuations at this scale we 
have that the result of this averaging is just the average 
density, i. e. one. Thus we obtain 



7i(0, 0) « exp - 



w[t',qit',x)] .(16) 



The above expression is the same expression that was 
derived in [15] for the singular density coarse-grained 
over the scale Idif- Using the result for the statistics 
of the coarse-grained density obtained in fl5|. we con- 
clude that the single-point statistics of the density is log- 
normal. 



JL) 

Idif J 



Cky pip-1) 



(17) 



We note that the fact that Idif is defined only up to 
a factor of order one is of no consequence. Any addi- 
tional factor C of order one would produce one when 
Idif is substituted by Cldif in the above equation, due 
to C'-^'^^ ~ 1. We also substituted I by 77 since by def- 
inition of I at this scale there are no fluctuations and 
{rj/l)'"'^^ ~ 1. Above it is assumed that p is not too 
large: very high-order moments are determined by the 
optimal fluctuations for which isotropic compression of 
the mass occurs, and compressibility is not small in any 
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sense [3,[i]. Finally, we note that since rj/ldif ^ Pr^/"^ 
then the above answer can be rewritten as 

{nP) « Pr^^^p(p"i)/2, (18) 

where we again use the fact that factors of order unity 
become irrelevant after being raised to a small power. 

Similar consideration of the different point correlation 
functions of the density reveals that these correlation 
functions are given by the correlation functions of the sin- 
gular density ris coarse-grained over the scale Idif- The 
latter correlation functions were computed in [15], where 
they were shown to obey the log-normal statistics. For 
example the expression for the pair correlation function 
is modified from Eq. ([5]) to 

(n(O)n(a;)) « f — ^) , x « r;. (19) 

Though the above formula is obtained by interpolation 
from X ^ Idif to a; 3> Idif, it gives an answer which is 
correct approximately. This is again because the correla- 
tion function changes insignificantly when x changes by 
the order of magnitude due to Cky 1- This finishes 
the consideration of the steady state statistics of the den- 
sity of particles moving according to Eq. ([6|). We note 
that our analysis that assumed that there are no signifi- 
cant fluctuations of density already at the scale / ^ is 
self-consistent as it is clear from the obtained formulas. 

We showed that it is possible to describe the statistics 
of the density of the particles in real turbulence taking 
into account the Brownian motion of the particles. The 
result is quite expectable. The statistics of the density is 
roughly the same as the statistics of the coarse-grained 
density of the particles that do not perform the Brown- 
ian motion, provided the scale of coarse-graining is cho- 
sen to be equal to the diffusive scale Idif- Still several 



non-trivial points must be stressed. First and foremost, 
the Brownian motion is always present for the real par- 
ticles and clear understanding of its effects is necessary. 
As an application let us derive the criterium for the exis- 
tence of significant fluctuations of the density. The latter 
exist only if (n^) is significantly larger than one, i. e. 
Cky In Pr ^ 1 ■ We observe that the dependence on the 
Prandtl number is only logarithmic. Thus since Cky is 
proportional to St^ times a power of the Reynolds num- 
ber, then the latter must be substantial to have signifi- 
cant fluctuations dX St <^1, cf. [5]. 

The second reason why the analysis is important is 
methodological. We saw that a different method was 
used to derive the statistics than in the previous works. 
Since the statistics of the distribution of particles with 
no Brownian motion is an idealization, the singular dis- 
tribution obtained in that limit must be considered as 
the K — )■ limit of the distribution obtained here. In 
other words, the physically sound way of thinking of the 
singular measure Us supported on the strange attractor 
of the particles in space is as the limit 



Us = lim nin), (20) 

where «(«;) is the smooth steady-state density obtained 
in this work, where k, is considered as a parameter. 

Our analysis can begeneralized directly to the case of 
light particles as in jl5| . A possible future application 
of the methods developed here is to the study of the 
distribution of inertial particles a,i St ^ 1, and it is the 
subject for the future work. 
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